Using multifrontal hierarchically solver and HPC systems for 3D Helmholtz problem

نویسندگان

  • Sergey Solovyev
  • Dmitry Vishnevsky
  • Hongwei Liu
چکیده

We present a multi-frontal hierarchically semi-separable solver to perform forward modeling of the 3D Helmholtz acoustic problem. Our frequency-domain solver combines two efficient approaches. First, it uses an optimal 27-point finite-difference scheme to decrease numerical dispersion and reduce required discretization of the model in terms of points per wavelength from 15 to about 4. Second, it uses a supernodal multi-frontal method based on low-rank approximation and hierarchically semi-separable (HSS) structure to improve performance, decrease memory usage and make it practical for realistic size 3D models required by full-waveform inversion. We performed validation and performance testing our new solver using a 3D synthetic model. Performance and OMP scalability of the solver were compared with the Intel MKL PARDISO.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Anisotropic ‘helmholtz’ Equations: Massively Parallel Structured Multifrontal Solver Using Nested Dissection Based Domain Decomposition with Separators of Variable Thickness

Abstract. We consider the discretization and approximate solution of inhomogeneous anisotropic ‘Helmholtz’ equations in 3D. The anisotropy comprises general (tilted) TI symmetries. In particular, we are concerned with solving these equations on a large domain, for a large number of different sources. We make use of a nested dissection based domain decomposition in a massively parallel multifron...

متن کامل

On 3D modeling of seismic wave propagation via a structured parallel multifrontal direct Helmholtz solver

We consider the modeling of (polarized) seismic wave propagation on a rectangular domain via the discretization and solution of the inhomogeneous Helmholtz equation in 3D, by exploiting a parallel multifrontal sparse direct solver equipped with Hierarchically Semi-Separable (HSS) structure to reduce the computational complexity and storage. In particular, we are concerned with solving this equa...

متن کامل

Efficient Scalable Algorithms for Hierarchically Semiseparable Matrices

Hierarchically semiseparable (HSS) matrix algorithms are emerging techniques in constructing the superfast direct solvers for both dense and sparse linear systems. Here, we develope a set of novel parallel algorithms for the key HSS operations that are used for solving large linear systems. These include the parallel rank-revealing QR factorization, the HSS constructions with hierarchical compr...

متن کامل

Superfast Multifrontal Method for Large Structured Linear Systems of Equations

In this paper we develop a fast direct solver for large discretized linear systems using the supernodal multifrontal method together with low-rank approximations. For linear systems arising from certain partial differential equations such as elliptic equations, during the Gaussian elimination of the matrices with proper ordering, the fill-in has a low-rank property: all off-diagonal blocks have...

متن کامل

Superfast Multifrontal Method for Structured Linear Systems of Equations

In this paper we develop a fast direct solver for discretized linear systems using the multifrontal method together with low-rank approximations. For linear systems arising from certain partial differential equations such as elliptic equations we discover that during the Gaussian elimination of the matrices with proper ordering, the fill-in has a low-rank property: all off-diagonal blocks have ...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2015